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Abstract 

We show that a twist of a three-dimensional tube of uniform cross- 
section yields an improved decay rate for the heat semigroup associated 
with the Dirichlet Laplacian in the tube. The proof employs Hardy in- 
equalities for the Dirichlet Laplacian in twisted tubes and the method of 
self-similar variables and weighted Sobolev spaces for the heat equation. 
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1 Introduction 



It has been shown recently in [7] that a local twist of a straight three-dimensional 
tube flo ■.— M.XLLioi non-circular cross-section w C leads to an effective repul- 
sive interaction in the Schrodinger equation of a quantum particle constrained 
to the twisted tube Qg . More precisely, there is a Hardy-type inequality for the 
particle Hamiltonian modelled by the Dirichlet Laplacian — A^" at its thresh- 
old energy Ei if, and only if, the tube is twisted (c/ Figure [1]). That is, the 
inequality 

-A'^'-E,>g (1.1) 

holds true, in the sense of quadratic forms in L'^{flg), with a positive function g 
provided that the tube is twisted, while g is necessarily zero for fig- Here Ei 
coincides with the first eigenvalue of the Dirichlet Laplacian — A"^ in the cross- 
section UJ. 





Figure 1 : Untwisted and twisted tubes of elliptical cross-section. 

The inequality (jl.ip has important consequences for conductance properties 
of quantum waveguides. It clearly implies the absence of bound states (i.e., 
stationary solutions to the Schrodinger equation) below the energy Ei even if 
the particle is subjected to a small attractive interaction, which can be either 
of potential or geometric origin (c/ [7] for more details). At the same time, a 
repulsive effect of twisting on eigenvalues embedded in the essential spectrum 
has been demonstrated in [14] . Hence, roughly speaking, the twist prevents 
the particle to be trapped in the waveguide. Additional spectral properties of 
twisted tubes have been studied in HSl [2] . 

It is natural to ask whether the repulsive effect of twisting demonstrated 
in [7] in the quantum context has its counterpart in other areas of physics, 
too. The present paper gives an affirmative answer to this question for systems 
modelled by the diffusion equation in the tube fig: 

Ut~Au^O, (1.2) 

subject to Dirichlet boundary conditions on dilg. Indeed, we show that the twist 
is responsible for a faster convergence of the solutions of (|1.2p to the (zero) stable 
equilibrium. The second objective of the paper is to give a new (simpler and 
more direct) proof of the Hardy inequality (jl.ip under weaker conditions than 
those in (Tj. 
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1.1 The main result 

Before stating the main result about the large time behaviour of the solutions 
to (|1.2p . let us make some comments on the subtleties arising with the study of 
the heat equation in Qg. 

The specific deformation Qg of fig via twisting we consider can be visu- 
alized as follows: instead of simply translating uj along R we also allow the 
(non-circular) cross-section to to rotate with respect to a (non-constant) angle 
xi 9{xi). See Figure [T] (the precise definition is postponed until Section [2 
cf Definition 1 2. 111 . We assume that the deformation is local, i.e., 



6 has compact support in K. (1-3) 

Then the straight and twisted tubes have the same spectrum {cf [T71 Sec. 4]): 

(7(-Ag'') =ae..(-Ag«) = [Si,oo). (1.4) 

The fine difference between twisted and untwisted tubes in the spectral setting 
is reflected in the existence of (|f .ip for the former. 

In view of the spectral mapping theorem, the indifference p.4p transfers to 
the following identity for the heat semigroup: 

W>0, ||e<^*|L.(^^)^,.(,,^)=e-^^*, (1.5) 

irrespectively whether the tube Qg is twisted or not. That is, we clearly have 
the exponential decay 

\Ht)\\mns) < e'""'' hoWmns) (1-6) 

for each time t > and any initial datum uq of p.2p . To obtain some finer dif- 
ferences as regards the time-decay of solutions, it is therefore natural to consider 
rather the "shifted" semigroup 

S{t) := e^^i+Ei)t (17) 

as an operator from a subspace of L'^{ng) to L^{fle). 

In this paper we mainly (but not exclusively) consider the subspace of initial 
data given by the weighted space 

L^{ng,K) with if(a;) := e^i/\ (1.8) 

and study the asymptotic properties of the semigroup via the decay rate defined 

by 

Tiftg) :=sup{r| 3Cr >0,Vi>0, ||^(i)IU=(f2„K)^L2(f2,) < Cr (1 + t)"^} . 

Our main result reads as follows: 

Theorem 1.1. Let 9 e C\R) satisfy ([O]). We have 



= 1/4 if fig is untwisted, 
> 3/4 if ^g is twisted. 
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The statement of the theorem for solutions u of (|1.2p in Qg can be reformu- 
lated as follows. For every F < F(r2g), there exists a positive constant Cr such 
that 

\\uit)\\mns) <Ct{1 + t)-^ e-^i* \\uo\\mne,K) (1-9) 

for each time t > and any initial datum uq E L^{il.g, K). This should be 
compared with the inequality (|1.6p which is sharp in the sense that it does 
not allow for any extra polynomial-type decay rate due to (jl.Sp . On the other 
hand, we see that the decay rate is at least three times better in a twisted tube 
provided that the initial data are restricted to the weighted space. 

A type of the estimate ()1.9p in an untwisted tube can be obtained in a less 
restrictive weighted space (c/ Theorem 14. ip . The power 1/4 actually reflects 
the quasi-one-dimensional nature of our model. Indeed, in the whole Euclidean 
space one has the well known dimensional bound 

W>0, ||e<*|L.(«.^^)^,.(«.) < (l+t)-'/'. (1.10) 

The fact that the power 1/4 is optimal for untwisted tubes can be established 
quite easily by a "separation of variables" (c/ Proposition [321) • The fine effect 
of twisting is then reflected in the positivity of F(rie) — 1/4; in view of p.lOp . 
it can be interpreted as "enlarging the dimension" of the tube. 



1.2 The idea of the proof 

The principal idea behind the main result of Theorem 1 1.1) i.e. the better decay 
rate in twisted tubes, is the positivity of the function g in (II. ip . In fact. Hardy 
inequalities have already been used as an essential tool to study the asymptotic 
behaviour of the heat equation in other situations f3','^. However, it should be 
stressed that Theorem 11.11 does not follow as a direct consequence of ()l.ip by 
some energy estimates (c/ Section H75|) but that important and further technical 
developments that we explain now are needed. Nevertheless, overall, the main 
result of the paper confirms that the Hardy inequalities end up enhancing the 
decay rate of solutions.. 

Let us now briefly describe our proof (as given in Section [5]) that there is 
the extra decay rate if the tube is twisted. 

I. First, we map the twisted tube fig to the straight one Qq by a change of 
variables, and consider rather the transformed (and shifted by Ei) equation 

ut-idi-ddryu-A'u-Eiu^O (1.11) 

in r^o instead of (|1.2p . Here —A' := — 5| — 9| and dr := X3d2 — X2d3, with 
X — {xi, X2, xs) € flo, denote the "transverse" Laplace and angular-derivative 
operators, respectively. 

II. The main ingredient in the subsequent analysis is the method of self-similar 
solutions developed in the whole Euclidean space by Escobedo and Kavian [5]. 
Writing 

u{yi,y2,y3, s) = e''/'^M(e^/^yi, ?/2, ys, - 1) , (1.12) 
the equation (jl.lip is transformed to 

Us-^yi diu - {di - Us drfu - l^'u- Eie" u- \ u^Q (1.13) 
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in self-similarity variables (y, s) G fio x (0, cxd), where 



(1.14) 



Note that p.l3p is a parabohc equation with time- dependent coefficients. 
This non-autonomous feature is a consequence of the non-trivial geometry we 
deal with and represents thus the main difficulty in our study. We note that an 
analogous difficulty has been encountered previously for a convection-diffusion 
equation in the whole space but with a variable diffusion cocfhcient 5 . 

III. We reconsider (jl.lSp in the weighted space (jl.Sp and show that the as- 
sociated generator has purely discrete spectrum then. Now a difference with 
respect to the self-similarity transformation in the whole Euclidean space is 
that the generator is not a symmetric operator if the tube is twisted. However, 
this is not a significant obstacle since only the real part of the corresponding 
quadratic form is relevant for subsequent energy estimates (c/ (|5.11|) ). 

IV. Finally, we look at the asymptotic behaviour of ()1.13p as the self-similar 
time s tends to infinity. Assume that the tube is twisted. The scaling coming 
from the self-similarity transformation is such that the function ()1.14p converges 
in a distributional sense to a multiple of the delta function supported at zero as 
s — > oo. The square of as becomes therefore extremely singular at the section 
{0} X w of the tube for large times. At the same time, the prefactors in (|1.13p 
diverge exactly as if the cross-section of the tube shrunk to zero s — > oo. Taking 
these two simultaneous limits into account, it is expectable that (jl.lSp will be 
approximated for large times by the essentially one-dimensional problem 

ips-\yiipy^-ipy^y^-\Lp = Q, s G ( , CX) ) , 1 £ R , (1.15) 

with an extra Dirichlet boundary condition at j;i = 0. This evolution equation 
is explicitly solvable in (R, K) and it is easy to see that 

\\v\\l'^{v.,k) < e"^'' W^oWl^r.k) , (1-16) 

for any initial datum ipQ. Here the exponential decay rate transfers to a poly- 
nomial one after returning to the original time t, and the number 3/4 gives rise 
to that of the bound of Theorem 1 1.1 1 in the twisted case. 

On the other hand, we get just 1/4 in (|1.16p provided that the tube is 
untwisted (which corresponds to imposing no extra condition at yi =0). 

Two comments are in order. First, we do not establish any theorem that 
solutions of (|1.13p can be approximated by those of (|1.15p as s ^ oo. We only 
show a strong-resolvent convergence for operators related to their generators 
(Proposition 15.4]) . This is, however, sufficient to prove Theorem 11.11 with help 
of energy estimates. Proposition 15.41 is probably the most significant auxiliary 
result of the paper and we believe it is interesting in its own right. 

Second, in the proof of Proposition 15.41 we essentially use the existence of 
the Hardy inequality (jl.ip in twisted tubes. In fact, the positivity of g is di- 
rectly responsible for the extra Dirichlet boundary condition of p.l5p . Since the 
Hardy inequality holds in the Hilbert space L^(r2o) (no weight). Proposition [53] 
is stated for operators transformed to it from ()1.8p by an obvious unitary trans- 
form. In particular, the asymptotic operator ho oi Proposition 15.41 acts in a 
different space, L^(M), but it is unitarily equivalent to the generator of (|1.15p . 
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1.3 The content of the paper 

The organization of this paper is as follows. 

In the following Section [2] we give a precise definition of twisted tubes 
and the corresponding Dirichlet Laplacian — A^®. 

Section [3] is mainly devoted to a new proof of the Hardy inequality (Theo- 
rem [311]) as announced in [T5]. We mention its consequences on the stability of 
the spectrum of the Laplacian fProposition 13 . 2| ) and emphasize that the Hardy 
weight cannot be made arbitrarily large by increasing the twisting (Proposi- 
tion 13. 3|) . Finally, we establish there a new Sobolev-type inequality in twisted 
tubes (Theorem [321) • 

The heat equation in twisted tubes is considered in Section [31 Using some 
energy-type estimates, we prove in Theorems 14.11 and 14.21 polynomial-type de- 
cay results for the heat semigroup as a consequence of the Sobolev and Hardy 
inequalities, respectively. Unfortunately, Theorem 14.21 does not represent any 
improvement upon the 1 /4-decay rate of Theorem l4 . 1 1 which is valid in untwisted 
tubes as well. 

The main body of the paper is therefore represented by Section [SI where we 
develop the method of self-similar solutions to get the improved decay rate of 
Theorem 11.11 as described above. Furthermore, in Section 15.91 we establish an 
alternative version of Theorem 11.11 

The paper is concluded in Section [51 by referring to physical interpretations 
of the result and to some open problems. 

2 Preliminaries 

In this section we introduce some basic definitions and notations we shall use 
throughout the paper. 

2.1 The geometry of a twisted tube 

Given a bounded open connected set w C M^, let fig := K x ai be a straight tube 
of cross-section uj. We assume no regularity hypotheses about uj. Let 6* : R ^ M 
be a C^-smooth function with bounded derivative (occasionally we will denote 
by the same symbol 9 the function 6* (g) 1 on SIq). We introduce another tube of 
the same cross-section w as the image 

where the mapping Cg ■.M.'^ ^ R"^ is given by 

Cg{x) :— (^xi,X2 cos6'(a;i) -I- sin6'(a;i), —X2 sin 6{xi) -f 2:3 cos6'(a;i)) . (2.1) 

Definition 2.1 (Twisted and untwisted tubes). We say that the tube VLg is 
twisted if the following two conditions are satisfied: 

1. 9 is not constant, 

2. w is not rotationally symmetric with respect to the origin in M^. 
Otherwise we say that Q.g is untwisted. 
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Here the precise meaning of uj being "rotationally symmetric with respect 
to the origin in M^" is that, for every ■& e (0, 27r), 

ujf) :— {x2 cos-d + xz sini^, —X2 sin-d + cos-d | (x2, X3) G cj} = tiJ , 

with the natural convention that we identify to and uj§ (and other open sets) 
provided that they differ on a set of zero capacity. Hence, modulus a set of zero 
capacity, lo is rotationally symmetric with respect to the origin in if, and 
only if, it is a disc or an annulus centered at the origin of M^. In view of this 
convention, any untwisted ile can be identified with the straight tube fio by an 
isometry of the Euclidean space. 

We write x = (xi, 2:2, 2:3) for a point/vector in . If x is used to denote 
a point in fip or fie, we refer to xi and x' := {x2,X3) as "longitudinal" and 
"transverse" variables in the tube, respectively. 

It is easy to check that the mapping Ce is injective and that its Jacobian 
is identically equal to 1. Consequently, Cg induces a (global) diffeomorphism 
between ilg and ^Ig. 

2.2 The Dirichlet Laplacian in a twisted tube 

It follows from the last result that Q.g is an open set. The corresponding Dirichlet 
Laplacian in L'^{VLg) can be therefore introduced in a standard way as the self- 
adjoint operator — A^" associated with the quadratic form 

Ql^m :=l|V*||i.(o,), *el)(gg«) -.^Hline). 

By the representation theorem, -A^"* = -A^" for 'J e ©(-A^") := {* e 
i?Q(r2e) I A^* S L^(r2e)}, where the Laplacian A^* should be understood in the 
distributional sense. 

Moreover, using the diffeomorphism induced by Cg, we can "untwist" the 
tube by expressing the Laplacian — A^" in the curvilinear coordinates deter- 
mined by (|2.ip . More precisely, let Ug be the unitary transformation from 
L^i^g) to L'^iVto) defined by 

C/e^t^fo^e. (2.2) 

It is easy to check that Hg :— Ug{—l\^)Ug^ is the self-adjoint operator in 
L^(r2o) associated with the quadratic form 

OeM \\di4'-0drnl^no) + \\^'nl^n,) > V' S ^{Qe) Hl{n^) . (2.3) 

Here V :— (82,83) denotes the transverse gradient and dr is a shorthand for 
the transverse angular-derivative operator 

dr -.^ T ■ V' = X3d2 - X2d3 , whcrc T(a;2, X3) := (xs, -0:2) . 

We have the point-wise estimate 

\drip\ < a |VV| , where a := sup \x'\ . (2.4) 

x' ^UJ 

The sesquilinear form associated with Qg[-] will be denoted by Qg{-, •). In the 
distributional sense, we can write 

Hgij^~{di-edr)^ij-A'^, (2.5) 

where —A' := — 9| — denotes the transverse Laplacian. 
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3 The Hardy and Sobolev inequalities 



In this section we summarize basic spectral results about the Laplacian — A^* we 
shall need later to study the asymptotic behaviour of the associated semigroup. 

3.1 The Poincare inequality 

Let El be the first eigenvalue of the Dirichlet Laplacian in uj. Using the Poincare- 
type inequality in the cross-section 

l|v/||i.(^)>£;i||/||i.(^), v/ei/iH, (3.1) 

and Fubini's theorem, it readily follows that Qei['0] ^ ^illV'lli2(s-;|,) for every 
ip G 7?Q(rio)- Or, equivalently, 

- Ag» > El (3.2) 

in the form sense in L^(r2e). Consequently, the spectrum of — A^" does not start 
below El . The result p.2p can be interpreted as a Poincare- type inequality and 
it holds for any tube ^Iq. 

The inequality p.2p is clearly sharp for an untwisted tube, since (|1.4p holds 
in that case trivially by separation of variables. In general, the spectrum of 
— A^* can start strictly above Ei if the twisting is effective at infinity (c/ [TSl 
Corol. 6.6]). In this paper, however, we focus on tubes for which the energy Ei 
coincides with the spectral threshold of — A^**. This is typically the case if the 
twisting vanishes at infinity (c/ [T71 Sec. 4]). More restrictively, we assume ()1.3p . 
Under this hypothesis, ()1.4p holds and (13. 2p is sharp in the twisted case too. 



3.2 The Poincare inequahty in a bounded tube 

For our further purposes, it is important that a better result than ()3.2p holds 
in bounded tubes. 

Given a bounded open interval / C M, let Hg be the "restriction" of Hg 
to the tube I x uj determined by the conditions di^p — 9drip = on the new 
boundary (dl) x oj. More precisely, Hg is introduced as the self-adjoint operator 
in L^{I X uj) associated with the quadratic form 

^e^iQi) -.^{^Pliixn) I ^ei/i(l]o)} . 

That is, we impose no additional boundary conditions in the form setting. 

Contrary to Hg, Hg is an operator with compact resolvent. Let X{6,I) 
denote the lowest eigenvalue of the shifted operator Hg — Ei. We have the 
following variational characterization: 

X{0,I)^ mm — ^ '-. (3.3) 

^eS(Q^)\{o} mhiixu) 

As in the unbounded case, (|3.ip yields that X{0,I) is non-negative (it is zero if 
the tube is untwisted). However, thanks to the compactness, now we have that 
Hg — ii'i is a positive operator whenever the tube is twisted. 



8 



Lemma 3.1. Let 6 G C^(R). Let / C M &e a bounded open interval such that 
6 \ L is not constant. Let uj he not rotationally invariant with respect to the 
origin in . Then 

X{0,I) > 0. 

Proof. We proceed by contradiction and assume that X{0,L) — 0. Then the 
minimum p. 311 is attained by a (smooth) function tp e D{Qg) satisfying (re- 
call dSl])) 

\\d^^P-9drij\\h^j^^)=0 and |lVV|li.(,,^)-£;i||^||i.(,,^)=0. (3.4) 

Writing tpi^) = 'Pixi)J^i{x') + 4'{x), where Ji is the positive eigenfunction 
corresponding to Ei of the Dirichlet Laplacian in L'^{uj) and (J7i, 0(a;i, ■))l'^(uj) — 
for every xi G /, we deduce from the second equality in (|3.4p that 4> = 0. The 
first identity is then equivalent to 

Ml^I)\\Jl\\lH..)HOv\\Ul)\\drM^^^^ 

Since (J7i, dT>Ji)L^(ui) = by an integration by parts, it follows that ip must be 
constant and that 

||0|U2(,)=O or |ia,Ji||L2(^) =0. 

However, this is impossible under the stated assumptions because ||0||2,2(/-| van- 
ishes if and only if 9 is constant on /, and drJ^i = identically in uj if and only 
if Lu is rotationally invariant with respect to the origin. □ 

Lemma 13.11 was the cornerstone of the method of [TH] to establish the exis- 
tence of Hardy inequalities in twisted tubes (see also the proof of Theorem 13.11 
below) . 

3.3 Infinitesimally thin tubes 

It is clear that A(^,R) := mia{He) = whenever (fOI holds {e.g., if (fOj) is 
satisfied) . It turns out that the shifted spectral threshold diminishes also in the 
opposite asymptotic regime, i.e. when the interval /e := (— e, e) shrinks, and this 
irrespectively of the properties of uj and 9. 

Proposition 3.1. Let 9 e C^{R). We have 

lim\(9,I,) =0. 

Proof. Let {i^k}kLo be an exhaustion sequence of oj, i.e., each ujk is an open set 
with smooth boundaries satisfying uJk (e cok+i and U^gWfc = cu. Let J^i denote 
the first eigenfunction of the Dirichlet Laplacian in L^(wfc); we extend it by zero 
to the whole R^. Finally, set := (1 (g) J^) o Ce^ with 9o{xi) ^(0) xi, i.e., 

''p'^ix) = Ji (^X2 cos(^oa;i) + 2:3 sin(^oa;i), -2:2 sin(^oa;i) + X3 cos(^oa;i)^ , 

where ^0 := ^(0). 
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For any (large) k E N there exists (smah) positive such that ijj'^ belongs 
to Tl{Qg') for all e < e^. Hence it is an admissible trial function for (|3.3p . Now, 
fix /c e N and assume that e < Cfc. Then we have 

U'\\lHi.><^)^\I^\\\Ji\\hM, 

where we have used the change of variables y = Cogix). At the same time, em- 
ploying consecutively the identity ditp'^ —6 drip'^ = (Oq—O) drip^ , the bound \2A\ . 
the identity \\/ ' -ip^ {x)\ = | V J'/^(?/2, J/a)! and the same change of variables, we get 
the estimate 

where the supremum norm clearly tends to zero as e — > 0. Finally, 

livvili.(,^x^) " ^^iliv>1li.(/.x.) = (^1 - E,) \h\ , 

where Ei denotes the first eigenvalue of the Dirichlet Laplacian in L'^{ujk)- 
Sending e to zero, the trial-function argument therefore yields 

lim A(^, /,) < - El . 

Since k can be made arbitrarily large and E'^ — > i?i as fc — > oo by standard 
approximation arguments (see, e.g., [1]), we conclude with the desired result. □ 

Remark 3.1 (An erratum to [17J). The study of the infinitesimally thin tubes 
played a crucial in the proof of Hardy inequalities given in 117] . According 
to Lemma 6.3 in [17j . A(0, /<:), with constant 6, is independent of e > (and 
therefore remains positive for a twisted tube even if e — > 0). However, in view 
Proposition 13.11 this is false. Consequently, Lemmata 6.3 and 6.5 and Theo- 
rem 6.6 in [17 cannot hold. The proof of Hardy inequalities presented in |17| 
is incorrect. A corrected version of the paper [17 can be found in [18] . 

3.4 The Hardy inequality 

Now we come back to unbounded tubes fig. Although p.2p represents a sharp 
Poincare-type inequality both for twisted and untwisted tubes (if (jl.4p holds), 
there is a fine difference in the spectral setting. Whenever the tube fig is non- 
trivially twisted (c/ Definition l2.ip . there exists a positive function g (necessarily 
vanishing at infinity if p.4p holds) such that (|3.2p is improved to (jl.ip . A variant 
of the Hardy inequality is represented by the following theorem: 

Theorem 3.1. Let G C^(M) and suppose that has compact support. Then 
for every ^' G HQ{ne) we have 

- Eiimhin,) > CHllpnUn,), (3.5) 

where p{x) :~ 1/^1 + x1 and ch is a non-negative constant depending on 
and uj. Moreover, ch is positive if, and only if, Qg is twisted. 
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Proof. It is clear that the left hand side of p.Sp is non-negative due to p.2p . The 
fact that = if the tube is untwisted follows from the more general result 
included in Proposition 13.21 2 below. We divide the proof of the converse fact 
(i.e. that twisting imphes ch > 0) into several steps. Recall the identification 
of * e L'^ing) with V := C/e^' e L^i^o) via ([221). 

1. Let us first assume that the interval / := (inf supp^, supsupp^) is symmetric 
with respect to the origin of M. 

2. The main ingredient in the proof is the following Hardy-type inequality for a 
Schrodinger operator in M x a; with a characteristic-function potential: 

llpV^lli^(Oo) < 16 -h (2 + 64/|/n (3.6) 

for every ip e HQ^flo). This inequality is a consequence of the classical one- 
dimensional Hardy inequality J^x^'^\ip{xi)\'^dxi < 4 J^\ip{xi)\'^dxi valid for 
any ip g 7Jo(M\{0}). Indeed, following jjj Sec. 3.3], let rj be the Lipschitz 
function on M defined by ri{xi) := 2|xi|/|/| for \xi\ < \I\/2 and 1 otherwise 
in R (we shall denote by the symbol the function 77 1 on M x cj). For any 
V' G C^(f^o), let us write = rjiP + (1 - r))^;, so that {r]i;){-,x') e i7o^(M\{0}) 
for every x' € iv. Then, employing Fubini's theorem, we can estimate as follows: 

llP^lli^(no) < 2 / xr'|(r7V')(a:)|'dx-f 2|l(l-77)^|li.(^^,) 
<8||ai(77V^)||i.(f,,,)+2||^||i.(,,^) 

< 16 \\vdiM\hin„) + 16 II (5ir])^||i.(o„) + 2 UWh^i-,^) 

< 16 l|5iVlli.(Oo) + (2 + 64/|/n \mh(iy.^) • 

By density, this result extends to all tp G iJg(r2o) = ^{Qe)- 

3. By Lemma ism we have 

g,[V']-i?i||V'lli2(o„) > QM-Ei\mh^i.,^^ > x{(),i)\mh^iy^-, (3.7) 

for every ijj E D{Qg). Here the first inequality employs the trivial fact that the 
restriction to / x w of a function from T){Qg) belongs to D{Qg). Under the 
stated hypotheses, we know from Lemma |3. II that X(0,I) is a positive number. 
4- At the same time, for every ip £ 'S{Qg), 



Oe[V']-i?iiiv^lli2(noj 



>e 1151^-1122(00) + 



1 - e 

2 



\V^P{x)\'-E,\i;{x)\Hdx 



>f||5i^lli2(o„)-Y^«'^^i 11^^111^(00) 

> e 1191^111^(00) - ll^lli==(/) IIV'lli^(7xc.) (3.8) 

for sufficiently small positive e. Here the first estimate is an elementary Cauchy- 
type inequality employing (|2.4p and valid for all e G (0, 1). The second inequality 
in (|3.8p follows from (|3.ip with help of Fubini's theorem provided that e is 
sufficiently small, namely if e < (l -I- a^||^||/,oo|-Rj 
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5. Interpolating between the bounds (|3.7p and l|3.8p . and using p.6p in the 
latter, we finally arrive at 

Qe[V]-i?i||V'lli.(o,)> llpV'lli.(ao) 



2 16 



for every e S)((5e). It is clear that the last line on the right hand side of this 
inequality can be made non-negative by choosing e sufficiently small. Such an e 
then determines the Hardy constant := e/32. 

6. The previous bound can be transferred to L'^{^lg) via (|2.2p . In general, if the 
centre of / is an arbitrary point G R, the obtained result is equivalent to 



where p^oix) := 1/^1 + {xi - x'()'^. This yields ([XS)) with 

/ • l + xj 
CH ■= Cjf mm — — fT-- , 

xies. 1 + {xi- x^y 

where the minimum is a positive constant depending on x^. □ 

The Hardy inequality of Theorem 13.11 was first established under addi- 
tional hypotheses. The present version is adopted from [H], where other variants 
of the inequality can be found, too. 



3.5 The spectral stability 

Theorem 13.11 provides certain stability properties of the spectrum for twisted 
tubes, while the untwisted case is always unstable, in the following sense: 

Proposition 3.2. Let V be the multiplication operator in L^{flg) by a bounded 
non-zero non-negative function v such that v(x) |a;i|~^ as \xi\ — > oo. Then 

1. if fie is twisted with £ C^(M) and has compact support, then there 
exists £o > such that for all e < Eq, 



2. if rig is untwisted then, for all e > 0, 

Mai-A^^" -eV) < Ei 



Proof. The first statement follows readily from one part of Theorem 13.11 To 
prove the second property (and therefore the other part of Theorem 13.11 stating 
that c^f = if the tube is untwisted), it is enough to consider the case — 
and construct a test function ip from _ffQ(r2o) such that 

Pom ||VV.||i.(f,,,) - i?i||Vlli.(no) - £ < 
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for all positive e. For every n > 1, we define 

ll^nix) -.^ ipn{xi)Ji{x2,X3,) , (3.9) 

where Ji is the positive eigenfunction corresponding to Ei of the Dirichlet 
Laplacian in the cross-section w, normalized to 1 in L'^{u>), and 

<^n{.xi) := exp ■ (3-10) 

In view of the separation of variables and the normalization of J7i , we have 

where Vi{xi) := \\v{xi,-y/^Ji\\^2(^^y By hypothesis, vi € L-^(R) and the in- 
tegral ||t'i||Li(E) is positive. Finally, an explicit calculation yields ||<i5n||L2(R) 
^-i/4_ gy ^YiQ dominated convergence theorem, we therefore have 

Po[ipn] > -e ||wi||l1(R) . 

n — >oo 

Consequently, taking n sufficiently large and e positive, we can make the form 
Poiipn] negative. □ 

Since the potential V is bounded and vanishes at infinity, it is easy to see 
that the essential spectrum is not changed, i.e., CTcssC—A^" — eV) — [i?i,oo), 
independently of the value of e and irrespectively of whether the tube is twisted 
or not. As a consequence of Proposition 13.21 we have that an arbitrarily small 
attractive potential —eV added to the shifted operator — A^" — Ei in the un- 
twisted tube would generate negative discrete eigenvalues, however, a certain 
critical value of e is needed in order to generate the negative spectrum in the 
twisted case. In the language of [3T], the operator — A^* — Ei is therefore 
subcritical (respectively critical) if Qe is twisted (respectively untwisted). 



3.6 An upper bound to the Hardy constant 

Now we come back to Theorem 13.11 and show that the Hardy weight at the right 
hand side of (|3.5p cannot be made arbitrarily large by increasing 9 or making 
the cross-section w more eccentric. 

Proposition 3.3. Let 9 G C^(R) and suppose that 9 has compact support. Then 

CH < 1/2, 

where ch is the constant of Theorem \3.1\ 

Proof. Recall the unitary equivalence of — A^" and Hg given by (|2.2p . We 
proceed by contradiction and show that the operator Hg ~ Ei ~ cp^ is not non- 
negative if c > 1/2, irrespectively of properties of 9 and uj. (Recall that p was 
initially introduced in Theorem 13.11 as a function on fig . In this proof, with an 
abuse of notation, we denote by the same symbol analogous functions on Qq 
and M.) It is enough to construct a test function ip from T){Qg) such that 

Pg'ii^] :=QeM-£;iiiV'iii.(^,„)-c||pviii2(^,„) <o. 
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As in the proof of Proposition 13. 21 we use the decomposition (|3.9p . but now the 
sequence of functions is defined as follows: 



V>n{xi) := < 



if xie[bi,bl). 
otherwise . 



Here {fe^jneN, with j = 1,2,3, are numerical sequences such that supsupp6' < 
< < ^or cach n G N and 6^ ^ cxd as n — > cxd; further requirements will 
be imposed later on. Since ipn and 9 have disjoint supports, and J7i is supposed 
to be normalized to 1 in L'^{uj), it easily follows that 

Note that the right hand side is independent of 6 and oj. An explicit calculation 
yields 

. 2 _ 1 1 

n n 

„ ,12 _ bl - bi + [[bif - l](arctan62 _ arctanfoi ) - bl log i±{^ 
IIP'^JL^W - (62 _ 51)2 ^ 

- ^« + [(^fJ^ - l](arctan&3 _ arctan62) _ bl log 



{bl-biy 



Specifying the numerical sequences in such a way that also the quotients b^/b\ 
and bl/b^i^ tend to infinity as n — > cxd, it is then straightforward to check that 

blPe^^n] ^l-2c. 

n — >oo 

Since the limit is negative for c > 1/2, it follows that Pg[ipn] can be made 
negative by choosing n sufficiently large. □ 

The proposition shows that the effect of twisting is limited in its nature, 
at least if (|1.3p holds. This will have important consequences for the usage of 
energy methods when studying the heat semigroup below. 

3.7 The Sobolev inequality 

Regardless of whether the tube is twisted or not, the operator — A^" — Ei 
satisfies the following Sobolev-type inequality. 

Theorem 3.2 (Sobolev inequality). Let 6 e Ci(R) and suppose that 9 has 
compact support. Then for every € -ffg (fig) H L'^(Q,e, p~^) we have 

\\'^nUn<.)-Ei\mhin,-> > ^^^lljp' (3.11) 



where := J^dx2dx3, (^J^dxi\{^ o Ce){x)\)'^ and cs is a positive constant 

depending on 9 and lo . 
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Proof. Recall that ^ o Ce = f/e* =: ip belongs to L'^i^o)- First of all, let 
us notice that ||5'||i is well defined for ^E" S {Q.g , p^^) . Indeed, the Schwarz 
inequality together with Fubini's theorem yields 

11*11? < IIP- Vlli.(oo) / = \\p-'nhin,) ^ < ^ • (3.12) 

Here the equality of the norms is obvious from the facts that the mapping Cg 
leaves invariant the first coordinate in and that its Jacobian is one. We also 
remark that, by density, it is enough to prove the theorem for ^ e C^{^g). 
The inequality (j3.1ip is a consequence of the one-dimensional inequality 

e H\R) n L\^) , > - , (3.13) 

^ IIv'IIli(r) 

which is established quite easily by combining elementary estimates 

< ||<y3||Li(R)ll<y5||L-(R) and ||¥'IIl~(R) <'^\\f\\L'^m\\^\\L^{WL)- 

In order to apply (|3.13p . we need to estimate the left hand side of p. lip from 
below by ||9i'0||^2(Q|j). We can proceed as in the proof of Theorem 13.11 Inter- 
polating between the bounds (|3.7p and (|3.8p . we get 

l|V*||i2(o,) - £^i||*lli.(o.) > I Wdinl^no) . 

where e =: 8 cs is a positive constant depending on 6 and lo. Using now p.l3p 
with help of Fubini's theorem, we conclude the proof with 

in ,||2 ^ 1 [ ||^(-,X2,a:3)||i.(K) 1 ll*lli.(n,) 



Here the second inequality follows by the Holder inequality with properly chosen 
conjugate exponents (recall also that |lV'||L2(n„) = ||*||l2(0£,))- D 

4 The energy estimates 
4.1 The heat equation 

Having the replacement u{x, t) i-^ e~^^* u{x, t) for p.2p in mind, let us consider 
the following t-time evolution problem in the tube fig: 

Ut — Au — Eiu = in ilg X (0, oo) , 

u = Mo in fie x {0} , (4.1) 
u — in ((9f2e) X (0, oo) , 

where uq G L^(rie). 

As usual, we consider the weak formulation of the problem, i.e., we say a 
Hilbert space-valued function u £ L'^^^(^{0, oo); Hq (fig)) , with the weak deriva- 
tive u' G L'^^^(^{0, oo); H^^ (Qg)) , is a (global) solution of (|4.ip provided that 

{v,u'it)) + {Vv,Vuit))^,^^^^~E, {v,uit))^,^^^^^ = (4.2) 
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for each v E ^[^{fle) and a.c. t E [0,oo), and u(0) = uo- Here (•, •) denotes the 
pairing of HQ{ilg) and H~^{ne). With an abuse of notation, we denote by the 
same symbol u both the function on f2e x (0, oo) and the mapping (0, cxd) — > 

Standard semigroup theory imphes that there indeed exists a unique solution 
of (|4.ip that belongs to C° (^[0 , oo) , L'^ (fig)) . More precisely, the solution is 
given by u{t) = S{t)uo, where S{t) is the heat semigroup (|1.7p associated with 
— A^" — El. By the Beurling-Deny criterion, S{t) is positivity-preserving for 
all t > 0. 

Since Ei corresponds to the threshold of the spectrum of — A^* if (|1.3p 
holds, we cannot expect a uniform decay of solutions of (|4.1[1 as i — s- oo in this 
case. More precisely, the spectral mapping theorem together with (|1.4p yields: 

Proposition 4.1. Let E C^(M) and suppose that has compact support. Then 
for each time t > we have 

\\S{t)\\L^(ae)-^L^(ng) = 1 • 

Consequently, for each t > and each s E (0, 1) we can find an initial datum 
uq E -ffo(^f) such that ||wo||l2(o^) = 1 and such that the solution of ()4.ip 
satisfies 

\\u{t)\\L^ne) > 1-e- 



4.2 The dimensional decay rate 

However, if we restrict ourselves to initial data decaying sufficiently fast at 
the infinity of the tube, it is possible to obtain a polynomial decay rate for 
the solutions of ()4.ip . In particular, we have the following result based on 
Theorem [ 



Theorem 4.1. Let E C^(M) and suppose that has compact support. Then 
for each time t > we have 



'5'(i)||L2(ne,p-2)_,L2(Os) < + — ^ 



where cs is the positive constant of Theorem \3.'2\ and p is introduced in Theo- 
remlSli 



Proof. The statement is equivalent to the following bound for the solution u 
of jn]): 

/ 4c, \"'/^ 

Vte[0,(X3), \\u{t)\\mn,) < \\p-^uo\\mne}h + ^t] , (4-3) 

where uq E L"^ {VIq , p^'^) is any non-trivial datum. It is easy to see that the 
real and imaginary parts of the solution of (|4.ip evolve separately. Further- 
more, since S{t) is positivity-preserving, given a non- negative datum uq, the 
solution u[t) remains non-negative for all i > 0. Consequently, establishing 
the bound for positive and negative parts of u(t) separately, it is enough to 
prove (|4.3p for non-negative (and non-trivial) initial data only. Without loss of 
generality, we therefore assume in the proof below that > for alH > 0. 
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Let {ipn}^^i be the family of moUifiers on M given by p.lOp : we denote by 
the same symbol the functions (pn (81 1 on M x D fig. Inserting the trial 
function 

Vn{x;t) := Lpri{xi)Un{x2,X'i]t) , Un{x2,Xs,]t) := 1 1 (^„u(- , 0:2 , ; 1 1 ^^j^^ , 

into ()4.2p . we arrive at (recall the definition of || • ||i from Theorem 13. 2p 

< (aif„(t),aiu(0)i2(o,) 

< \\dlVn{t)\\L^(^ne)\\'^u{t)\\L2(n,) . 



Here the first inequality is due to the Poincare-type inequality in the cross- 
section p.ip . We clearly have 

\\dlVn{t)\\L^no) = ll<y3n||L2(K) 1 1 (i) 1 1 L2 = WVnWmR) 1 1 <^n"(i) 1 1 1 ■ 

Integrating the differential inequality, we therefore get 



/ \\^u{t')\\hin,)dt' . 

Since |l</5n||L2(R) and {Vn}'^=i is an increasing sequence of functions con- 
verging pointwise to 1 as n ^ oo, we conclude from this inequality that 

Vte[0,oo), < lluolli, (4.4) 

where ||wo||i is finite due to (|3.12p . 

Now, substituting u for the trial function v in Ij4.2p . applying Theorem 
and using (|4.4p . we get 



< -cs ■ 



WuoWt ■ 

An integration of this differential inequality leads to 



-1/4 



Dividing the last inequality by \\p ^ua\\L2(ne) ^-nd replacing ||wo||i with 
||p"^Mo||L2(na) using (|3.12p, we get 



»H^<,(l,l^«.,)-"'<(l,l^A-'\ 

\\P 'wo||L2(Oe) V 7^ / V ^ / 

where ^ ||wo||L2(ij(,)/||/0"^uo||L2(n(,) e (0, 1). This establishes (g^l). □ 
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As a direct consequence of the theorem, we get: 

Corollary 4.1. Under the hypotheses of Theorem \4-l\ T{^e) > 1/4. 

Proof. It is enough to reahze that L^{flg, K) is embedded in L'^{Q.g, P'^)- 

The foUowing proposition shows that the decay rate of Theorem 14. II is opti- 
mal for untwisted tubes. 

Proposition 4.2. Let U,g be untwisted. Then for each time t > we have 

\\SmLHn„K)-,mne) > ^ (1 + t)"'/' . 

Proof. Without loss of generality we may assume 9 = 0. It is enough to find 
an initial datum uq E L? (fio , K) such that the solution u of (|4.ip satisfies 

Vte[0,oo), > ^A^^^r'" ■ (4.5) 

The idea is to take uq := ipm where {i/'nlJ^Li the sequence p.9p approximating 
a generalized eigenfunction of — corresponding to the threshold energy Ei. 
Using the fact that is a cross-product of M and oj, (14. ip can be solved explicitly 
in terms of an expansion into the eigenbasis of the Dirichlet Laplacian in the 
cross-section and a partial Fourier transform in the longitudinal variable. In 
particular, for our initial data we get 



Htr 



where the second equality is a result of an explicit calculation enabled due to 
the special form of (pn given by (|3.10p . At the same time, for every n < 8 
tpn belongs to L^{ilo,K) and an explicit calculation yields 



12 ^2 



For the special choice n = 6 we get that the left hand side of (|4.5p actually 
equals the right hand side with t being replaced hy 2t/3 < t. □ 



The power 1 /4 in the decay bounds of Theorem 14.11 and Proposition 
reflects the quasi-one-dimensional nature of ilg [cf (jl.lOp ). at least if the tube 
is untwisted. More precisely. Proposition 14. 21 readilv implies that the inequality 
of Corollary 2T] is sharp for untwisted tubes. 

Corollary 4.2. Let VLg be untwisted. Then T{Q.g) = 1/4. 

This result establishes one part of Theorem 11.11 The much more difficult 
part is to show that the decay rate is improved whenever the tube is twisted. 
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4.3 The failure of the energy method 

As a consequence of combination of direct energy arguments with Theorem l3.11 
we get the foUowing resuh. In Remark 14.21 below we explain why it is useless. 

Theorem 4.2. Let ftg be twisted with 9 G C^(K). Suppose that 6 has compact 
support. Then for each time t > we have 

ll^(t)llL^(n«,p-)^L^(a.) < (i + 2t)-"^'"{i/^-«/^> , (4.6) 
where ch is the positive constant of Theorem \3.1\ 

Proof. For any positive integer n and x £ ile, let us set p„(a;) := mm{p{x), n^^}. 
Then is a non-decreasing sequence of bounded functions converging 

pointwise to p"""^ as n — > oo. Recalling the definition of p from Theorem 13.11 
it is clear that x t-^ Pn{x) is in fact independent of the transverse variables x' . 
Moreover, p~''u belongs to H}^{Vlo) for every 7 e M provided u € H}^{^e)- 

Choosing v := Pn'^u in ()4.2p (and possibly combining with the conjugate 
version of the equation if we allow non-real initial data), we get the identity 

\jt\\Pn'<tW = -\\p-'Vu{t)f+E,\\p;Mt)\?-'^{u{t)Vp-\Vu{t)) . (4.7) 

Here and in the rest of the proof, || -H and (•, •) denote the norm and inner product 
in L^(f2e) (we suppress the subscripts in this proof). Since p„ depends on the 
first variable only, we clearly have V(p^^) = (— 2/9^'^9ip, 0, 0). Introducing an 
auxiliary function Vn{t) :— p~^u{t), one finds 

\\p-^Vu{t)\\' = ||Vw„(i)||' + \\idip.n/pn)Vn{t)f +'2fi{yn{t),{diPn/pn)diVn{t)^ , 
?Si(uit)Vp-^,Vu{t)^ =^2\\{diPn/pr,)Vn{t)f-^'ii{vn{t),{dlPn/Pn)diVn{t)^ . 

Combining these two identities and substituting the explicit expression for p, 
we see that the right hand side of (|4.7p equals 

-\\S/Vn{t)\\^ + Ei\\Vn{tW + \\{diPn/pn)Vn{tW 

- -||Vt;„(t)f + E,\\Vn{t)\\^ + WXnPVnitW - WXnP^Vnim^ (4.8) 
< (1 - ch) \\XnPVn{t)f - ||XnP'«n(Oll' ■ 

Here Xn denotes the characteristic function of the set ilg :— ilg n {supp(9ip„)}, 
and the inequality follows from Theorem l3.1l and an obvious inclusion fig C fig. 
Substituting back the solution u{t), we finally arrive at 

lj^\\Pn'u{t)\\^ < (1 - Ch) WXnPVnim^ ~ II XnP'^^n (0 f 

<{^-CH)\\XnPVn{t)f ■ (4.9) 

Now, using the monotone convergence theorem and recalling the initial data 
to which we restrict in the hypotheses of the theorem, the last estimate implies 
that u{t) belongs to L^{flg,p~^) and that it remains true after passing to the 
limit n oo, i.e., 

~\\p-Mm'<ii-CH)Mt)f. (4.10) 
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At the same time, we have 



~\Ht)r = -{\\^<t)r-EiMt)r) 



2 



< -CH \\pu{t)\ 

<-c.^^, (4.11) 
\\P ^u{t)P 

where the equahty follows from (j4.ip . the first inequality follows from Theo- 
rem 13.11 and the last inequality is established by means of the Schwarz inequal- 
ity- 
Summing up, in view of (|4.1ip and (|4.10p . a{t) and b{t) := 
\\p^^u{t)\\'^ satisfy the system of differential inequalities 

a<-2cH — , 6<2(l-cff)a, (4.12) 



with the initial conditions a(0) = ||uo|P ='■ oo and 5(0) ~ ||p^^uo|P =: bo. We 
distinguish two cases: 

1. ch >l- In this case, it follows from the second inequality of (|4.12p that b 
is decreasing. Solving the first inequality of (|4.12p with b being replaced by bg, 
we then get 

a{t) <ao[l + 2cH {aa/bo)t]~\ 

Dividing this inequality by bo and maximizing the resulting right hand side with 
respect to ao/6o G (0, 1), we finally get 

Vte[0,oo), \\u{t)\\ < \\p-\o\\{l + 2cHt)-^^^ , (4.13) 

which in particular implies (|4.6p . 

2. ch < I- We "linearize" (|4.12p by replacing one a of the square on the right 
hand side of first inequality by employing the second inequality of (|4.12p : 

a „ a ch b 

<-2cHT< 



a b 1 — Ch b 

This leads to 

a/ao < {b/boY^^ . 

Using this estimate in the original, non-linearized system (|4.12p . i.e. solving the 
system by eliminating b from the first and a from the second inequality of (|4.12p , 
we respectively obtain 

a{t) < ao [1 + 2 {ao/bo) t] , b{t) < 5o [l + 2 (ao/6o) <] ■ 

Dividing the first inequality by bo and maximizing the resulting right hand side 
with respect to ao/bo G (0, 1), we finally get 

Vie[0,oo), \\uit)\\ < \\p-\oUl + 2t)-'"/' , (4.14) 
which is equivalent to (|4.6p . □ 
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Remark 4.1. We see that the power m the polynomial decay rate of Theorem 
diminishes as ch — *■ 0. Let us now argue that this cannot be improved by the 
present method of proof. Indeed, the first inequality of (j4.1ip is an application 
of the Hardy inequality of Theorem 13.11 and the second one is sharp. The 
Hardy inequality is also applied in the first inequality of (|4.9|) . In the second 
inequality of (|4.9p . however, we have neglected a negative term. Applying the 
second inequality of (|4.1ip to it instead, we conclude with an improved system 
of differential inequalities 

2 2 

a<-2ci/^, &< 2(l-Cij)a-2^. (4.15) 



The corresponding system of differential equations has the explicit solution 
00(^7777 ^TT— ^1 ' KO = «W(l + VF[eoexp(^o + 2t)]), 

where ^0 ^0/^0 ^ 1 > and W denotes the Lambert W hmction (product 
log), i.e. the inverse function of w 1— > ?i;exp(u'). Since 

W [Co exp (Co + 2t)] =2t + o{t) as t ^ 00 , 

we see that the t^"^"/^ decay in (|4.6p for < 1 cannot be improved by replac- 
ing (|iT^ with gTH). 

Remark 4.2. Note that the hypothesis (|1.3p is not explicitly used in the proof of 
Theorem l4.21 it is just required that the inequality p.Sp holds with some positive 
constant ch- For tubes satisfying (jl.Sp . however, we know from Proposition [33] 
that the constant cannot exceed the value 1/2. Consequently, irrespectively of 
the strength of twisting, Theorem 14.21 never represents an improvement upon 
Theorem 14. II This is what we mean by the failure of a direct energy argument 
based on the Hardy inequality of Theorem 13.11 



5 The self-similarity transformation 

Let us now turn to a completely different approach which leads to an improved 
decay rate regardless of the smallness of twisting. 



5.1 Straightening of the tube 

First of all, we reconsider the heat equation (|4.ip in an untwisted tube fip by 
using the change of variables defined by the mapping Cg . In view of the unitary 
transform (|2.2p . one can identify the Dirichlet Laplacian in L^(ile) with the 
operator (|2.5p in L^(rio), and it is readily seen that (|4.ip is equivalent to 

ut + Hgu — Eiu = in ^Iq x (0, 00) , 

plus the Dirichlet boundary conditions on dflo and an initial condition at i = 0. 
(We keep the same latter u for the solutions transformed to ilo.) More precisely, 
the weak formulation (|4.2p is equivalent to 

{v,u'it)) + Qg{v,uit)) ~ E,{v,u{t)) (5.1) 
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for each v G HQ{no) and a.e. t E [0,oo), with u(0) ^ uq E L'^{ilo). Here (•, •) 
denotes the pahmg of HQ{no) and i7~^(S7o)- We know that the transformed 
solution u belongs to C°([0, oo), L^(ilo)) by the semigroup theory. 

5.2 Changing the time 

The main idea is to adapt the method of self-similar solutions used in the case of 
the heat equation in the whole Euclidean space by Escobedo and Kavian [8J to 
the present problem. We perform the self-similarity transformation in the first 
(longitudinal) space variable only, while keeping the other (transverse) space 
variables unchanged. 

More precisely, we consider a unitary transformation U on L^(rio) which 
associates to every solution u E ifoc((Oj oo), dt; L^(rio, dx)) of (|5.ip a self-similar 
solution u := Uu in a new s-time weighted space ifoc((Oi '^)^ c^ds; L^{Vl[), dy)) 
via ()1.12p . The inverse change of variables is given by 

u{xi,X2,X^, t) = (t+ l)-l/4 u{{t + iy'^/^Xi,X2,X3,\0g{t + 1)) . 

When evolution is posed in that context, y = 2/2, 2/3) plays the role of space 
variable and s is the new time. One can check that, in the new variables, the 
evolution is governed by ()1.13|) . 

More precisely, the weak formulation (j5.1|) transfers to 

{v,u'{s) - ^yidiu{s)) +Qs{i,u{s)) - Ei {v,u{s)) ^^^^^^^^ =0 (5.2) 

for each v E HQ{flo) and a.e. s E [0, 00), with m(0) — uq := Uuo = uq. 
Here Qsi', •) denotes the sesquilinear form associated with 

Qs[u] := \\diu - as 9rw||i2(Oo) + || V'u||i2(f2„) - ^ ||M||i2(Oo) , 

uEDiQs) ■.= H^{no), 

where as has been introduced in (|1.14p . 

Note that the operator Hs in L'^i^o) associated with the form Qs has s-time- 
dependent coefficients, which makes the problem different from the whole-space 
case. In particular, the twisting represented by the function p.l4p becomes 
more and more "localized" in a neighbourhood of the origin yi = for large 
time s. 

5.3 The natural weighted space 

Since U acts as a unitary transformation on L^(r2o), it preserves the space norm 
of solutions of (|5.ip and (|5.2p . i.e., 

\\u{t)\\L^no) = Ms)\\L^no) ■ (5.3) 

This means that we can analyse the asymptotic time behaviour of the former 
by studying the latter. 

However, the natural space to study the evolution (|5.2p is not i^(rio) but 
rather the weighted space (jl.Sp . For fc e Z, we define 

Hk -.^L^ {no, K\y,)dyidy2dy3). 
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Hereafter we abuse the notation a bit by denoting by initiaUy introduced as 
a function on fig in (jl.Sp . the analogous function on R too. Note that K^^^^ is 
the first eigenfunction of the harmonic-osciUator Hamiltonian 

/^:=-^ + ^2/? in L\R) (5.4) 

[i.e. the Friedrichs extension of this operator initiaUy defined on C^(R)). The 
advantage of reformulating (|5.2p in Tii instead of Ho = L^(rio) lies in the fact 
that then the governing elliptic operator has compact resolvent, as we shall see 
below (c/ Proposition [531) • 

Let us also introduce the weighted Sobolev space 

nl := Hl{n„, K\yi) dyidy2dys) , 

defined as the closure of C^{flo) with respect to the norm (|1 ■ llwt + 11^ ' II Wt)^^^- 
Finally, we denote by Ti,^^ the dual space to Tij.. 



5.4 The evolution in the weighted space 

We want to reconsider ()1.13p as a parabolic problem posed in the weighted 
space Til instead of Ho- We begin with a formal calculation. Choosing v{y) :— 
K{yi)v{y) for the test function in (|5.2p . where v G C^{flo) is arbitrary, we can 
formally cast (|5.2p into the form 

{v,u'{s)) +as{v,ii{s)) ^0. (5.5) 

Here (•, •) denotes the pairing of Til and Ti-i^, and 

as{v, u) :— (^div — ag drV, diu — CTs dru)^^ + (V'w, V'u)^^ 

- El {v, u)^^ - i {yi v, a, dru)^^ - ^ {v, u)^^ . 

Note that Ug is not a symmetric form. 

Of course, the formulae are meaningless in general, because the solution u{s) 
and its derivative u'{s) may not belong to Hi and Hi^, respectively. We there- 
fore proceed conversely by showing that (|5.5p is actually well posed in Hi and 
that the solution solves (|5.2p too. As for the former, we have: 

Proposition 5.1. For any uq e Hi, there exists a unique function u such that 

u e LL((0,(X3);7ii) n C"([0,oo);Hi) , u' e ^^((0, «.); Tif') , 

and it satisfies (|5.5p for each v G Hi and a.e. s € [0,oo), and u(0) = uq. 

Proof. First of all, let us show that is well-defined as a sesquilinear form 
with domain D{as) '■= Hi for any fixed s G [0,oo). In view of the boundedness 
of (Ts (for every finite s) and the estimate (|2.4p . it only requires to check that 
yivGHi provided v G Hi. Let v G C^{flo). Then 

\\yivrni=2f yiHy)\'^^dy 

= -2 f {\viy)\^ + 2yi^[viy)diviy)]}Kiyi)dy 
< 4||2/if||wi ||<9iw|l-Hi • 
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Consequently, 

\\yiv\\m<M\div\\ni<M\v\\^i. (5.6) 

By density, this inequality extends to all v &Hl. Hence, as{v, u) is well defined 
for all s > and v,u € Hi (we suppress the tilde over u in the rest of the 
proof). Then the Proposition follows by a theorem of J. L. Lions [H Thm. X.9] 
about weak solutions of parabolic equations with time-dependent coefficients. 
We only need to verify its hypotheses: 

1. Measurability. The function s i-^ as{v,u) is clearly measurable on [0, oo) for 
all v,u £ T-L\, since it is in fact continuous. 

2. Boundedness. Let sq be an arbitrary positive number. Using the bound- 
edness of 6, the estimates ()2.4p and (|5.6p . it is quite easy to show that there 
is a constant C, depending uniquely on sq, H^lliooj-R) and the geometry of uj 
(through a and Ei), such that 

\as{v,u)\<C\\v\\ni\\u\\ni (5.7) 

for all ,s G [0, sq] and v, u E Hi- 

3. Coercivity. Recall that Us is not symmetric and that we consider complex 
functional spaces. However, since the real and imaginary parts of the solution u 
of (|5.5p evolve independently, one may restrict to real-valued functions v and u 
there. Alternatively, it is enough to check the coercivity of the real part of a^. 
We therefore need to show that there are positive constants e and C such that 
the inequality 

^{as[v]}>e\\v\\l,i-C\\v\\i^^ (5.8) 
holds for all v £ Til and s £ [0, sq], where as[v] := as{v, v). We have 

3?{a,H} = \\d,v - a. drv^n^ + e/ ||V'H|?,, - E, \\vf^^ - ^ \\v\\l^^ 

- ^di{yiv,asdrv)Hi (5.9) 
for all V £ Hi- For every v £ C^{ilQ), an integration by parts shows that: 

n{yiv,asdrv)ni^O; (5.10) 

by density, this result extends to all v £ Ti\. Hence, the mixed term in ()5.9p 
vanishes. We continue with estimating the first term on the right hand side 
of (15:911: 



\div - as c'tw||«i > e - ^ Iks drv\\l^^ 



valid for every e £ (0, 1) and v £ Hi- Here the second inequality follows from 
the definition of cr^ in (|1.14p and the estimate (|2.4p . Using p.ip with help of 
Fubini's theorem, wc therefore have 
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provided that e is sufficiently small (so that the expression in the round brackets 
is positive). Putting this inequality into (|5.9p . recalling (|5.10p and using the 
trivial bounds 1 < e'* < e*" for s S [0, sq], we conclude with 

^{as[v]} > e||Vi;| 

valid for all sufhciently small e and all real- valued v e Ti-l. It is clear that the 
last inequality can be cast into the form (|5.8p , with a constant e depending on a 
and ||0||ioo(u), and a constant C depending on sq, ||0||ioo(K) and the geometry 
of w (through a and Ei). 

Now it follows from pT, Thm. X.9] that the unique solution u of (|5.5p satisfies 

i2e i2((0,so);H}) nC°([0,so];Hi) , u' G L2((0,so);Hr') ■ 

Since sq is an arbitrary positive number here, we actually get a global continuous 
solution in the sense that u S C''([0,oo);7ii). □ 

Remark 5.1. As a consequence of (|5.7p . (|5.8p and the Lax-Milgram theorem, it 
follows that the form as is closed on its domain Hi- 

Now we are in a position to prove a partial equivalence of evolutions (|5.2p 
and ([53]) . 

Proposition 5.2. Lei uq G Tii. Lei w &e the unique solution to (j5.5p /or eac/i 
w G 7i} and a.e. s G [0,cx)), subject to the initial condition u(0) — uq, that is 
specified in Provosition \5. 11 Then u is also the unique solution to (j5.2p for each 
V G Tip and a.e. s G [0, cx)), subject to the same initial condition. 

Proof. Choosing v{y) :— K{yi)^^ v{y) for the test function in (|5.5p . where v G 
C|5"(rio) is arbitrary, one easily checks that u satisfies (|5.2p for each v G C^(rio) 
and a.e. s G [0,cx3). By density, this result extends to all w G Hq. □ 




5.5 Reduction to a spectral problem 

As a consequence of the previous subsection, reducing the space of initial data, 
we can focus on the asymptotic time behaviour of the solutions to (|5.5p . Choos- 
ing V := u{s) in ()5.5p (and possibly combining with the conjugate version of the 
equation if we allow non-real initial data), we arrive at the identity 

~\m\\ni^-jPH^)], (5.11) 

where ^{as[u]}, u G D(ji^^) := ©(a^) = Til (independent of s). 

Recalling (|5.9p and (|5.10p . we have 

Jfi) [u] = \\diu - a, drufn, + e^ || V'^iH^,^ - i^i e^ \\u\\l,^ - J ll^iH?,^ . 

As a consequence of (|5.7p . (15.81) and the Lax-Milgram theorem, we know that 
ji^'' is closed on its domain 7i}. It remains to analyse the coercivity of the 
form Js^\ 
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More precisely, as usual for energy estimates, we replace the right hand side 
of (|5.1ip by the spectral bound, valid for each fixed s e [0,oo), 

yneHl, Ji'\u]>fi{s)\\uf^^, (5.12) 

where fi{s) denotes the lowest point in the spectrum of the self-adjoint opera- 
tor Ti^^ in Hi associated with Js^\ Then (|5.1ip together with ()5.12p implies 
the exponential bound 

Vs e [0, oo) , \\u{s)\\n, < WuoWn, e" ^(r)dr ^ ^r, -^^^ 

In this way, the problem is reduced to a spectral analysis of the family of oper- 
ators {ri^''}5>o. 

5.6 Removing the weight 

In order to investigate the operator Ts^^ in Tii, we first map it into a unitarily 
equivalent operator Ti°^ in Hq. This can be carried out via the unitary transform 
Uo -Ti.! Ho defined by 

{Uou){y) ■.^K^/\y,)u{y). 

We define rj"' UqTs^'^U^^ ^ which is the self-adjoint operator associated with 
the quadratic form Ji"^ [u] := ji^\u^^v], v € :^ UoSi{//^). A straight- 

forward calculation yields 

- \\div - a.dMl^i, + ^ Wvivfno + W^'vWho - Eie^ • (5.14) 

It is easy to verify that the domain of Js^^ coincides with the closure of 
C^ino) with respect to the norm (|| • \\^^^ + \\V ■ \\^^ + \\yi ■ In par- 

ticular, ©(jj^-*) is independent of s. Moreover, since this closure is compactly 
embedded in Ho (one can employ the well-known fact that (|5.4p has purely 
discrete spectrum, which essentially uses the fact that the form domain of h is 
compactly embedded in L^(K)), it follows that Ti°^ (and therefore Ts^^) is an 
operator with compact resolvent. In particular, we have: 

Proposition 5.3. Tg^^ ~ Tg^^ have purely discrete spectrum for all s £ [0, oo). 
Consequently, fi{s) is the lowest eigenvalue of Ts^\ 

5.7 The asymptotic behaviour of the spectrum 

In order to study the decay rate via (|5.13p , we need information about the limit 
of the eigenvalue fi(s) as the time s tends to infinity. 

Since the function as from p.l4p converges in the distributional sense to 
a multiple of the delta function supported at zero as s oo, it is expectable 
(c/ (|5.14p ) that the operator ri^-* will converge, in a suitable sense, to the one- 
dimensional operator h from (|5.4p with an extra Dirichlet boundary condition 
at zero. More precisely, the limiting operator, denoted by ho, is introduced as 
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the self-adjoint operator in L^(R) whose quadratic form acts in the same way 
as that of h but has a smaUer domain 

^{h]/') := e V{h^/') I ^(0) = 0} . 

1 /2 

Ahernatively, the form domain 'S{hj^ ) is the closure of C^(R \ {0}) with 
respect to the norm (|| • Wl^^^^ + ||V • + \\yi ■ HisfR))^/^. 

To make this limit rigorous (Tj"' and ho act in different spaces), we fol- 
low [ini and decompose the Hilbert space Ha into an orthogonal sum 

where the subspace Sji consists of functions of the form t/jiiu) = f{yi)J'i{y')- 
Recall that J7i denotes the positive eigenfunction of — A'J^ corresponding to Ei, 
normalized to 1 in L'^{uj)- Given any ip £ Hq, we have the decomposition 
= V'l + with -01 G -^1 as above and (p £ ^j^. The mapping tt : ip i-^ i/j i is 
an isomorphism of L^(]R) onto ^i. Hence, with an abuse of notations, we may 
identify any operator h on i^(R) with the operator tt/itt"^ acting on i^i C Hq. 

Proposition 5.4. Let be twisted with 9 e C^(R). Suppose that 9 has com- 
pact support. Then Tg^^ converges to hjj 0"*^ in the strong-resolvent sense as 
s ^ CO, i.e., for every F € Ho, 



lim 



{T^"^ + l)-'F-\{hD + iy'(BO^ 



F 



■Ho 

± 



Here denotes the zero operator on the subspace S^^ C Ho- 

Proof. For any fixed F G Ho and sufficiently large positive number z, let us set 
Tps ■= {Ts^^ -f z)^^F. In other words, satisfies the resolvent equation 

Vwe2)(Ji"'), 4''Hv,ibs) + ziv,ibs)Ho = iv,F)no. (5.15) 
In particular, choosing ^pg for the test function v in (|5.15p . we have 

= {i^s,F)no<\Us\\'no + \\F\\no- (5-16) 

Henceforth we assume that z > 1/4. 

We employ the decomposition ips{y) = fs{yi)Ji{yi) + 4'siy) where (ps £ Sjj^, 
i.e., 

VyieR, (j'i,(/.,(2/i,-))i.(^) -0. (5.17) 
Then, for every e € (0, 1), 

I! VV.|1^„ - E,U.,r^^ = e\\V<Psfn„ + (1 - e)|lV>.||?,„ - E^U^f^^ 
>e\\V<Ps\\'no + [i^-^-)E2-E,]\\<Ps\\ii„, 

where E2 denotes the second eigenvalue of — A'J^. Since Ei is (strictly) less 
then E2, we can choose the e so small that (j5.16p implies 

Us\\'n„<Ce-' and |1VV.||?,„ < ^^6"% (5.18) 
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where C is a constant depending on lo and Hi^Hwo- At the same time, (|5.16p 
yields 

||<Ps|1l2(r) < C, |lyi'^s|lL2(R) < C, and bi^sHwo < C* , (5-19) 

where C is a constant depending on 

To get an estimate on the longitudinal derivative of 'i/'s, we handle the first 
three terms on left hand side of (I5.16|) as follows. Defining a new function Us £ 
Tio by ipsiy) = e^^^Us{e^^^yi,y') {cf the self-similarity transformation p.l2p ) 
and making the change of variables {xi, x') = (e^^^yi, y'), we have 

Ji°)[^.] - e^\\d,u, -ed.usfn^^ + ^— Wx.UsWl,^^ + (|| V'^.,||2,„ - i^i 
> CH WpusWn^ , 

= e' chWpsiI'sWh^. where p,(y) := ^(e'/^yi, y') ■ (5-20) 

In the second inequality we have employed the Hardy inequality of Theorem l3.H 
the constant ch is positive by the hypothesis. Consequently, ()5.16p yields 

Msrn,<Ce-^, (5.21) 

where C is a constant depending on 6, uj and Now, proceeding as in the 

proof of p.Sp . we get 

for every e < (l + a^||^|||<^(jg))"\ where 1^ := e""/'^! = {e'^/'^xi \ xi£ 1} with 
/ := (inf supp sup supp 0). Since 

\\i^s\\L-(I^y.u.)<C\\psiJs\\H,, (5.22) 

where C is a constant depending exclusively on /, (|5.16p together with (|5.2ip 
implies H^iV'sII^q < C, where C is a constant depending on ^, cj and 
Recalling ()5.17p . we therefore get the separate bounds 

\\diC^s\\H,<C and ||<^.|U2(r) < C , (5.23) 

with the same constant C. 

By (|5.18p . (ps converges strongly to zero in TYq as s ^ oo. Moreover, it follows 
from dnHH]), (jgl^ and ([g?^ that {(1)s}s>q is a bounded family in I)(ji"^). 
Consequently, (ps converges weakly to zero in S)(ji*''') as s ^ oo. 

At the same time, it follows from (|5.19p and (|5.23p that {^Ps}s>o is a bounded 
family in T){h}/'^). Therefore it is precompact in the weak topology of X)(/i^/^). 
Let Lpoo be a weak limit point, i.e., for an increasing sequence of positive numbers 
{s„}„gN such that s„ ^ oo as n cxi, {ips^}n&i converges weakly to tpoo in 
I)(fti/2). Actually, we may assume that it converges strongly in (M) because 
S)(/i^/-^) is compactly embedded in _L^(R). 
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Employing (|5T7)) . fOT]) together with ([O^ gives 

ysWh^j^-, < Ce-\ (5.24) 

where C is a constant depending on 0, w and Multiplying this inequality 

by e*/^ and taking the limit s — s- cx), we verify that 

(^oc(0) = 0. (5.25) 

(We note that Dih^/"^) C H^{R) and that -ff^(J) is compactly embedded in 
C°'^{J) for every A G (0, 1/2) and any bounded interval J C K.) 

Finally, let e C(f (M\{0}) be arbitrary Taking v{x) := ¥3(xi)Ji(x') as the 
test function in (|5.15|) . with s being replaced by s„, and sending n to infinity, 
we easily check that 

'i'oo)L2(R) + — {yiip, yi^Poo)L^(M) + Z if, foo)L^{M) = if, /)l2(r) , 

where f{xi) := (Ji, -F(a;i, •))l2(w)- That is, ipao = (ho + z^^f, for any weak 
limit point of {fs}s>o- 

Summing up, we have shown that tps converges strongly to ip 
s->oo, where Voo(2/) := 9500(^1)^1(2/') = [{ho + z)-^ ® 0^]F. □ 

Remark 5.2. The crucial step in the proof is certainly the usage of the Hardy 
inequality in the second inequality of (|5.20p . Indeed, it enables one to control 
the mixed terms coming from the first term on the left hand side of (|5.16p . We 
would like to mention that instead of the Hardy inequality itself we could have 
used in (|5.20p the corner-stone Lemma [XT] This would leave to the lower bound 

ji°^[V's] > A(^, /) ||V's|li2(/^xw)' which is sufficient to conclude the proof in the 
same way as above. 

Corollary 5.1. Let fig be twisted with 9 e C^(R). Suppose that 9 has compact 
support. Then 

lim ^l{s) = 3/4. 



Proof. In general, the strong-resolvent convergence of Proposition 15.41 is not 
enough to guarantee the convergence of spectra. However, in our case, since the 
spectra are purely discrete, the eigenprojections converge even in norm (c/ 23J). 
In particular, ^(s) converges to the first eigenvalue of Hd. It remains to notice 
that the first eigenvalue of Hd coincides (in view of the symmetry) with the 
second eigenvalue of h which is 3/4. (For the spectrum of h, see any textbook 
dealing with quantum harmonic oscillator, e.g., [131 Sec. 2.3].) □ 

5.8 The improved decay rate - Proof of Theorem 11.11 

Now we have all the prerequisites to prove Theorcm ll.il Recall that the iden- 
tity r(r26() = 1/4 for untwisted tubes is already established by Corollary] 
Throughout this subsection we therefore assume that fig is twisted with 
and show that there is an extra decay rate. 

We come back to (|5.13p . It follows from Corollary 15.11 that for arbitrarily 
small positive number e there exists a (large) positive time such that for all 
s > Se, we have /i(s) > 3/4 — e. Hence, fixing e > 0, for all s > s^, we have 



/o 



^(r) dr <- fi{r) dr-{3/4 - £)(s - s^) < (3/4 - e)s£-(3/4 - e)s , 
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where the second inequality is due to the fact that /i(s) is non-negative for aU 
s > (it is in fact greater than 1/4, cf Proposition 15. 5p . At the same time, 
assuming e < 3/4, we triviaUy have 

- / n{r) dr<0< (3/4 - e)se-(3/4 - £)s 



also for all s < s^. Summing up, (j5.13p implies 

Ms)\\n,<C,e-^'/*-'>\\no\\n, (5.26) 

for every s € [0,oo), where := e*' > e'^^^'*"'^^''^ . Returning to the variables 
in the straightened tube via u = U'-^u, using (|5.3p together with the point- wise 
estimate 1 < K, and recalling that uq = uq, it follows that 

ll"Wllwo - Ms)\\no < \m\\n, <CAl + t)-^'/'-'^ hoik, 

for every t € [0,00). Finally, we recall that the weight K in Hi depends on the 
longitudinal variable only, which is therefore left invariant by the mapping Cg. 
Consequently, we apply the unitary transform (j2.2p and conclude with 



\\sit)U.^n,,K)^mn,) = sup <cAi + ty^'/^-^^ 

«oe-Hi\{o} \\uq\\hi 

for every t e [0,cx)). Since e can be made arbitrarily small, this bound implies 
r(f^0) > 3/4 and concludes thus the proof of Theorem ll.il 



5.9 The improved decay rate - an alternative statement 

Theorem 11.11 provides quite precise information about the extra polynomial 
decay of solutions u of p.2p in a twisted tube in the sense that the decay 
rate r(rie) is at least three times better than in the untwisted case. On the 
other hand, we have no control over the constant Cr in (jl.9p (in principle it 
may blow up as F — > F(Oe)). As an alternative result, we therefore present also 
the following theorem, where we get rid of the constant Cr but the prize we pay 
is just a qualitative knowledge about the decay rate. 

Theorem 5.1. Let 9 e C^{R) satisfy dO]). We have 

Vi>0, \\S{t)\\L2^n,M)^mn,) < {l + tr'^'''''^''\ (5.27) 

where ^ is a non-negative constant depending on 6 anduj. Moreover, 7 is positive 
if, and only if, is twisted. 

In order to establish Theorem 15.11 the asymptotic result of Corollary 15.11 
need to be supplied with information about values of /i(s) for finite times s. 



5.9.1 Singling the dimensional decay rate out 

It follows from Theorem 14.11 that there is at least a 1/4 polynomial decay rate 
for the solutions of the heat equations. In the setting of self-similar solutions 
(recall (j5.13p and the relation between the initial and self-similar times t and s 
given by (|1.12p ). this will be reflected in that we actually have /i(s) > 1/4, 
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regardless whether the tube is twisted or not. It is therefore natural to study 
rather the shifted operator Tj"'' — 1/4. However, it is not obvious from ()5.14p 
that such an operator is non-negative. 

In order to introduce the shift explicitly into the structure of the oper- 
ator, we therefore introduce another unitarily equivalent operator := 
U-iTs^\U-i)~^ in 7i_i, where the map U-i : Hq — > 7Y-i acts in the same 
way as Uq : 

{U^iv){y) ■.= K''^{y,)v{y). 

Ts^'^^ is the self-adjoint operator associated with the quadratic form ji^^'' [w] :— 
jf\{U^i)-'^'w\, w e £>(ji~^^) := Again, it is straightforward to 

check that 

Ji^i) [w] = \\d,w - a, drw\\i,_^ + e/ WVwfn^, E, + \ . 

Now it readily follows from the structure of the quadratic form that the 
shifted operator Tg — 1/4 is non- negative. Moreover, it is positive if, and 
only if, the tube is twisted. 

Proposition 5.5. If VIq is twisted with 9 €E C"'^(IR), then we have 

Vse[0,oo), /i(s)>l/4. 

Conversely, /i(s) = 1/4 /or all s G [0, oo) if fig is untwisted. 

Proof. Since ji ^\w] — j > for every w G we clearly have 

fJ-is) ^ 1/4, regardless whether the tube is twisted or not. By definition, if it 
is untwisted, then either a-g — identically in R for all s e [0,oo) or 9-rJi = 
identically in w, where J7i is the positive eigenfunction corresponding to Ei 
of the Dirichlet Laplacian in L^{tj). Consequently, choosing w{y) = Ji{y') as 
a test function for Js~^\ we also get the opposite bound /i(s) < 1/4 in the 
untwisted case. To get the converse result, we can proceed exactly as in the 
proof of Lemma [3. II Assuming /i(s) = 1/4 in the twisted case, the variational 
definition of the eigenvalue /i(s) would imply 

lks||L2(a,x-i) = or |19rJi||L2(^) = 0, 

a contradiction. □ 

Now we are in a position to prove Theorem 15. II 

5.9.2 Proof of Theorem [5A] 

Assume p.3p . It follows from Proposition l5.5l and Corollary El] that the number 

7:= inf u(s)-l/4 (5.28) 

se[o,oo) 

is positive if, and only if, is twisted. In any case, (|5.13|) implies 

\\u{s)\\n, < WuoWn, e-^^+^Z^)^ 
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for every s € [0,oo). Using this estimate instead of (|5.26p . but following the 
same type of arguments as in Section [HH] below (|5.26p . we get 



for every t € [0, oo). This is equivalent to (|5.27p and we know that 7 is positive 
if ^Iff is twisted. On the other hand, in view of Proposition [121 estimate (j5.27p 
cannot hold with positive 7 if the tube is untwisted. This concludes the proof 
of Theorem O 



6 Conclusions 

The classical interpretation of the heat equation (|1.2p is that its solution u 
gives the evolution of the temperature distribution of a medium in the tube 
cooled down to zero on the boundary. It also represents the simplest version 
of the stochastic Fokker-Planck equation describing the Brownian motion in Qg 
with killing boundary conditions. Then the results of the present paper can 
be interpreted as that the twisting implies a faster cool-down/death of the 
medium/Brownian particle in the tube. Many other diffusive processes in nature 
are governed by (|1.2p . 

Our proof that there is an extra decay rate for solutions of (|1.2p if the tube is 
twisted was far from being straightforward. This is a bit surprising because the 
result is quite expectable from the physical interpretation, if one notices that 
the twist (locally) enlarges the boundary of the tube, while it (locally) keeps 
the volume unchanged. (By "locally" we mean that it is the case for bounded 
tubes, otherwise both the quantities are infinite of course.) At the same time, 
the Hardy inequality (jl.ip did not play a direct role in the proof of Theorems ll.il 
and 15.11 (although, combining any of the theorems with Theorem 13. 11 we even- 
tually know that the existence of the Hardy inequality is equivalent to the extra 
decay rate for the heat semigroup). It would be desirable to find a more direct 
proof of Theorem 1 1 . 1 1 based on p.ip . 

We conjecture that the inequality of Theorem I 1 . 1 l ean be replaced by equality, 
i.e., r(rie) = 3/4 if the tube is twisted and ()1.3p holds. The study of the 
quantitative dependence of the constant 7 from Theorem [ST] on properties of 9 
and the geometry of lu also constitutes an interesting open problem. Note that 
the two quantities are related by 7 + 1/4 < r(ri6)). 

Throughout the paper we assumed p.3p . We expect that this hypothesis 
can be replaced by a mere vanishing of 9 at infinity to get Theorems 11.11 and 15.11 
(and also Theorem l3.ip . This less restrictive assumption is known to be enough 
to ensure (|1.4p and there exist versions of (jl.ip even if (|1.3p is violated (c/ [18]). 
However, it is quite possible that a slower decay of 9 at infinity will make the 
effect of twisting stronger. In particular, can T{Ug) be strictly greater than 3/4 
if the tube is twisted and 9 decays to zero very slowly at infinity? 

Equally, it is not clear whether Proposition 13.31 holds if (|1.3p is violated. 
There are some further open problems related to the Hardy inequality of The- 
orem [23] In particular, it is frustrating that the proof of the theorem does not 
extend to all 9 merely vanishing at infinity. In this context, it would be highly 
desirable to establish a more quantitative version of Lemma 13.1) i. e. to get a 
positive lower bound to X{9, 1) depending explicitly on 9, \I\ and oj. 
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On the other hand, a completely different situation will appear if one allows 
twisted tubes for which 9 does not vanish at infinity. Then the spectrum of 
— A^* can actually start strictly above Ei {cf [U] or [TTl Corol. 6.6]) and an 
extra exponential decay rate for our semigroup S{t) follows at once already 
in L'^(fle). In such situations it is more natural to study the decay of the 
semigroup associated with — A^" shifted by the lowest point in its spectrum. As 
a particularly interesting situation we mention the case of periodically twisted 
tubes, for which a systematic analysis based on the Floquet-Bloch decomposition 
could be developed in the spirit of [3D] . 

We expect that the extra decay rate will be induced also in other twisted 
models for which Hardy inequalities have been established recently [161 [E] ■ 

It would be also interesting to study the effect of twisting in other physical 
models. As one possible direction of this research, let us mention the question 
of the long time behaviour of the solutions to the dissipative wave equation 

Let us conclude the paper by a general conjecture. We expect that there 
is always an improvement of the decay rate for the heat semigroup if a Hardy 
inequality holds: 

Conjecture. Let ft be an open connected subset o/M'*. Let H and iJ+ be two 
self-adjoint operators in L^(0) such that 'mia{H) — infcr(iJ+) = 0. Assume 
that there is a positive smooth function : f2 — > R such that > g, while 
H ^ V is a negative operator for any non-negative non-trivial V G C^{n). 
Then there exists a positive function K : fl --*M. such that 

\\e'"+*\\L^n,K)^LHn) ^ ^ 
t^oo ||e--f^*||i2(s^^^)^i2(n) 

A similar conjecture can be stated for the same type of operators in different 
Hilbert spaces. In this paper we proved the conjecture for the special situation 
where H — — Ei and iJ+ = Hg — Ei (transformed Dirichlct Laplacians) in 
L^(ri), with = ilo (unbounded tube). In general, the proof seems to be a 
hardly accessible problem. 
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